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We present recent measurements by SLD of the branching fractions of B 
hadrons to states with and 2 open charm hadrons, BRqu and BR2D, from 
which both the average charm yield per B decay, N c and the inclusive branch- 
\ ing ratio into rare modes not containing any charmed hadrons, BR rare can be 

derived. We also present a new measurement of the B& mixing frequency Am^ 
and limits on the B s mixing frequency Am s . These analyses take advantage 
of the excellent vertexing resolution of the VXD3, a pixel-based CCD vertex 
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c3 detector, which enables the topological separation of the B and cascade D decay 

j^j ■ vertices. 

(N : 

C^ ■ 1 B Decay Charm Counting 

This measurement is motivated by the possible discrepancy between measurements and Standard 
Model predictions of the B semileptonic branching fraction and the average charm quark/antiquark 
yield per B decay, N c . Recent next-to-leading-order calculationsB agree with experiment only at 
rather low renormalization scales where it is not obvious that the perturbation series is nearing 
convergence. See Yamamotd for a discussion. A true discrepancy between the measurements 
and the Standard Model predictions could signal either a breakdown of standard calculational 
techniques, or an enhancement of the rate of decays into rare or unexpected new modes. 

<j_) ■ In this paper, we present measurements using a novel vertexing technique, of the branching 

ratios of B hadron decays into final states with or 2 weakly decaying charmed hadrons. These 
decay categories will be referred to as OD decays and 2D decays in which D refers to any of 
' D° / / D® /c-b&ryon. The sample of B hadrons used is the B^ /B° /B® /b— baryon admixture 

produced in decays of the Z°. The measurements of BRqd and BR2D may be combined with 
previous measurements of the B decay branching ratio into final states including charmonium, 

H ' BR c£ to get a experimental value for N c using: 

JV C = lx BR 1D + 2 x BR 2D + 2 x BR c5 

= 1.0 - BR QD + BR 2D + 2 x BR cS . (1) 

The second equation here is obtained by using BRqd + BRiu + BR 2 d = 1-0. The branching 
ratio into rare final states not containing any charm hadrons may also be calculated as: 

BRrare = BRqd — BR c5 (2) 

These measurements therefore not only yield the average charm count N c , but also elucidate the 
composition of the non-semileptonic decay width. In addition, the systematic uncertainties in 
these measurements are mostly uncorrelated with those of most previous measurements which 
have relied on counting D's via reconstruction of exclusive decay modes. These measurements 
therefore provide important new data on inclusive properties of B decays. 



"Invited talk presented at the XXXVIIth Rencontres de Moriond, Les Arcs, March 9-16, 2002. 
b Work supported in part by Department of Energy contract DE-AC03-76SF00515. 
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1.1 The Method 



The measurements of BRqd and BR2D described here utilize a correspondence between the 
number of heavy hadron weak decays in an event and the number of distinct topological decay 
vertices that can be reconstructed in the detector. A OD decay should produce a single secondary 
B decay vertex in addition to the primary Z° decay vertex. A ID decay should produce 
an additional secondary vertex at the D decay position, and a 2D decay should produce two 
additional secondary vertices. 

The correspondence between the number of heavy hadron weak decays and the number of 
reconstructed vertices is not perfect due to vertex finding inefficiencies due to low vertex track 
multiplicities and finite tracking resolution. Tails in the tracking resolution distribution may 
also cause false vertices to be formed away from the true decay positions, resulting in vertex 
finding overefficiencies. Due to these difficulties, it is not possible to identify the topological 
category of decays on an event-by-event basis. However, a counting analysis is still possible, 
using a suitable unfolding procedure to account on average for the effects of these vertexing 
issues. 

The nD branching ratios may be simultaneously measured by simply counting the number 
of vertices found in each B decay, and fitting the secondary vertex count (N v t x ) distribution 
for the entire data set to a linear combination of a set of distributions predicted by the MC for 
each decay category. In order for this procedure to work, both the MC detector and physics 
simulations must be carefully tuned to accurately model the vertex finding efficiency so that the 
correct probability distribution functions for N vtx for each decay category are produced. The 
MC detector simulation is calibrated using a number of supplementary measurements. The MC 
physics simulation is tuned to measurements made by MARKIII, CLEO, and LEP experiments. 

1.2 Fitting the vertexing distributions 

To form the vertex count distributions, first a sample of hadronic Z° decays is selected by 
requiring at least 7 measured tracks in the drift chamber and a visible energy of at least 30 
GeV in the calorimeters. The interaction point (IP) where the Z° decays is measured using 
the reconstructed tracks. Z° — > bb events are then selected by requiring at least one i?-taggecu 
hemisphere in each event. This tag yields a 98% pure sample of bb events. To get a sample of 
generic B decays, only the decays in the hemisphere opposite a^-B-tagged hemisphere are chosen. 
In the resulting sample of B decay hemispheres, the ZVTOPH ghost track algorithm is used to 
reconstruct the B and cascade D vertices. 

The vertex count distribution measured in the data is fit to a linear combination of distri- 
bution shapes predicted by the MC for each of the udsc background in the -B-tagged sample, 
the OD decays, the ID decays, and the 2D decays. The vertex count distribution shapes for 
each of these categories are shown in figure |l[ These shapes show significant differences which 
provide high analyzing power for the fit. The OD shape is strongly peaked in the 1-vertex bin 
as expected. Both a B and a D vertex are found about 50% of the time in the ID category, and 
a second D is also found a good fraction of the time in the 2D category. In this last case, due 
to low vertex track multiplicities, even an infinitely precise detector could find all three vertices 
at most ~ 40% of the time. 

In order to utilize extra available information from the measured vertex positions, these 
vertex count distributions are expanded as follows. The measured B decay length, defined as 
the distance between the measured IP position and the nearest reconstructed vertex in the 
hemisphere, is histogrammed separately for each value of the secondary vertex count. These 
histograms, shown in figure ^, can be viewed as slices of a 2-dimensional histogram with N v t x 
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on one axis and the measured B decay length on the other axis. These new shape distributions 
still provide information on the number of D's through N v t x and but now lifetime information 
from the long distance behavior of the decay length distributions as well as vertexing resolution 
information from the short distance behavior of the decay length distributions is also included. 
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Figure 1: MC predicted found secondary vertex Figure 2: MC predicted IP-vertex 1 separations (cm) 

count. Each histogram has been normalized to unit on a log scale. The rows represent JV S „=1,2,3,>4, and 

weight. the columns represent udsc background, 0D, ID, and 

2D. The histograms have been normalized to a data- 
sized sample. 



As before, the 2-dimensional vertex count distribution is fit to a linear combination of the 
2-dimensional distribution shapes predicted by the MC for each of the four decay categories. 
The fit is a binned x 2 fit- Each of the N v t x =1, 2, 3,> 4 decay length distributions is divided 
into 14 bins. In order to provide a balance between the information from short distance scales 
and that from long distance scales, the bins are chosen to be uniform in the logarithm of decay 
length. One extra bin is used for the N v t x = count where there is no explicit decay length 
measurement. The fitting function used is: 

where F £>, Fid, F2D, Fbkgd are the four normalized MC distributions, and i = {1..57} is the bin 
number. The parameters extracted from the fit are normalization R n , the udsc background 
fraction Rbkgd in the B tag and the branching ratios BRqd, BR2D- BR\e> has been eliminated 
to impose the constraint that the branching fractions sum to unity. 

The result of the fit, including the contributions to the measured shapes from the various 
sources, is shown in figure ||. The measured distribution appears to be modelled fairly well 
although the fit x 2 /d.o.f. = 1.6 is rather large. Since the detector resolution has been calibrated 
in several ways, the remaining discrepancies are believed to be due to imprecise modelling of the 
momentum spectrum of daughter particles at each decay stage. Variations of the MC modelling 
of the B decays are included in the systematic errors. 
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1.3 The Results 



The results of the measurement are: 



BR(B^(0D)X) = (3.7 ±1.1 ±2.1)% (4) 
BR{B^(2D)X) = (17.9 ± 1.4 ±3.3)% (5) 

where the first error is statistical and the second is systematic!. The correlation coefficients 
between the two measurements are Cqd,2D = 0.7^2 and —0.080 for statistical and systematic 



(6) 



errors, respectively. iV c is calculated using a valufl of BR^ = (2.3 ± 0.3)% in equation |T|: 



1.188 ± 0.010 ± 0.040 ± 0.006. 



Here, the third error is due to the uncertainty in BR cS . The measured value of N c is plotted in 
figure |H and compared with the LEP and CLEO measurement averages! (updated with the new 
estimates for and charmonium production!), and with the theory predicted region!. The 
plot indicates that the region of consistency is still at a low renormalization scale 

Limits on BR rare may be set using equation yielding a value of BR rare = (1.4 ± 2.4)% 
consistent with the theoretical expectation] of (2.6 ± 1.1) 
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Figure 3: Fit results on a log scale. The stacked MC 
histograms are from bottom to top: udsc background, 
0D, ID, and 2D. (a)-(d) show the IP to vertex 1 
separation [cm], (f)-(h) shows the nearest neighbor 
vertex separations [cm], and (e) shows the resulting 
match in the vertex count distribution. 



Figure 4: The SLD measurement of N c compared 
with measurements of N c and BRsl by LEP and 
CLEO and with theoretical expectations. The new 
SLD N c result is shown as a horizontal band on this 
plot. 



2 B Mixing 

The magnitude of the CKM matrix element Vtd may be extracted from measurements of the 
mixing frequencies Am^ and Am s . The ingredients of the mixing analyses include 1) a selection 
of neutral B decays from the -B-tagged sample described above; 2) an initial state flavor tag and 
3) a final state flavor tag in order to determine whether a B has mixed before decaying; and 4) a 
decay length and boost measurement to yield the proper time of the decay. For Am s , the mixing 
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Table 1: Comparison of Am s analyses in the number of events assigned to each analysis, the core and tail decay 
length resolution, the B s fraction in each event sample, and the final state mistag rate w. 





events 


ct l Om) 
core (tail) 


B s fraction 


w 


Charge Dipole 


11462 


81 (297) 


16% 


22% 


Lepton + D 


2087 


54 (213) 


16% 


4% 


Df + Tracks 


361 


50 (151) 


38% 


10% 



frequencies are then measured by performing amplitude fits on the plots of the mixed fraction 
as a function of proper time to extract the amplitude of each Fourier mode. A likelihood fit 
to the mixed fraction plot is used to extract Am^. Here we report the results of four analyses, 
three for Am s and one for Am^. More information may be found elsewheriH i. 

All four analyses share the same techniques for obtaining the neutral B sample, for doing 
the initial state tag, and for measuring the B boost. The charge of the decaying B is measured 
by associating charged tracks with a reconstructed B decay vertex to compute the net vertex 
charge. Selecting the zero charge bin yields a sample which is 87% true neutral decays. 

The initial state flavor tag is naturally provided by forward-backward asymmetry in the 
decay due to the 73% polarization of the SLC e~ beam. Other opposite hemisphere information 
such as the vertex charge, the momentum weighted jet charge, the charges of identified kaons 
and/or high P t leptons, and the dipole charge (see below), are also used in a neural network 
optimized tag to achieve a mistag rate of only 22%. 

The B boost is measured using a combination of calorimetry and of the kinematics of the 
reconstructed B vertex. a p /p « 0.08 for the 60% Gaussian core resolution, and 0.20 for the 
tail resolution. 

The techniques for the final state tag and decay length measurement are described below. 
For Am s the analyses are prioritized in the order below from highest to lowest sensitivity. Each 
B decay is assigned to the highest sensitivty analysis that it qualifies for. 



2.1 The Am s Analyses 

This l Dg + Tracks' method utilizes a fully reconstructed decay to identify the sign of the b 
quark. Requiring a Djr in the event enhances the B s fraction by rejecting many of the B^ decays. 
The reconstructed Djr trajectory is intersected with the remaining tracks from the B decay to 
get a precise decay length measurement. Because of resulting precise proper time resolution, 
this method contributes the most analyzing power at large Am s despite having low statistics. 

The 'Lepton+D' method uses a high-i^ identified lepton from the B decay to very cleanly 
identify the sign of the b quark. Selecting semileptonic decays rejects the decays with a wrong- 
sign D^r which dilute the tag purity of the other two analyses. The lepton trajectory is intersected 
with the inferred trajectory of a reconstructed D vertex in the event to measure the B decay 
length. 

The most inclusive analysis uses the 'Charge Dipole' technique which exploits the naturally 
occurring charge dipole in the B s — > DjX + decay cascade. In hemispheres with both the B and 
the D vertex reconstructed, the charge dipole may be defined as the charge difference (Qd — Qb) 
between the two vertices multiplied by the vertex separation distance. B s (B s ) decays will tend 
to have negative (positive) values of this quantity. The proper time is then measured from the 
decay length, defined as the distance between the IP and the closest secondary vertex. 
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The results of these three analyses are combined to form the amplitude fit plot shown in 
figure ||. Based on this plot, a lower bound of Am s < ll.lps -1 at 95% confidence level is 
calculated. This bound gives an upper bound on \Vtd\ in the p — rj plane of the CKM matrix. 

2.2 The Am,j Analysis 

This analysis tags the sign of the b quark using the identified from the cascade D decay. This 
tag is calibrated with the data by simultaneously fitting for Am^ and the 'right sign fraction' in 
the B decay sample (figure ||). The latter fraction is measured to be (79.7 ±2.2)% of the sample 
of 7844 B decays used. The mixing frequency is measured to be Arrid = 0.503 ± 0.028 ± 0.020. 
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Figure 5: Ami likelihood fit to the mixed 
fraction vs. proper time. 
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Figure 6: Combined amplitude fit results for Am s . 
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